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The impurity problems within vortex cores of two-dimensional s-wave and chiral p-wave su- 
perconductors are studied numerically in the framework of the quasiclassical theory of super- 
conductivity and self-consistent Born approximation under a trial form of the pair potential. 
The dispersion and impurity scattering rate (the inverse of the relaxation time) of the Andreev 
bound state localized in vortex cores are deduced from the angular-resoloved local density of 
states. The energy dependence of the impurity scattering rates depends on the pairing symme- 
try; particularly, in the chiral p-wave vortex core where chirality and vorticity have opposite sign 
and hence the total angular momentum is zero, the impurities are ineffective and the scattering 
rate is vanishingly small. Owing to the cancellation of angular momentum between chirality and 
vorticity, the chiral p-wave vortex core is similar to locally realized s-wave region and therefore 
non-magnetic impurity is harmless as a consequence of Anderson's theorem. The results of the 
present study confirm the previous results of analytical study (J. Phys. Soc. Jpn. 69 (2000) 
3378) in the Born limit. 

KEYWORDS: superconductivity, vortex, quasiclassical theory, impurity effect, self-consistent Born approximation, 
chiral superconductor, flux flow conductivity 
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§1. Introduction 

Vortex cores are often regarded as normal regions 
with the radius of the coherence length £ - This pic- 
ture is valid only in the dirty superconductors where 
the mean free path I is much smaller than £n. In clean 
superconductors where I 3> £o, quasiparticles experi- 
ence the Andreev reflections 1 -' much often than colli- 
sions with impurities. The constructive interference at 
the multiple Andreev reflections leads to the Andreev 
bound states 2 - 1 in the cores; these states are nothing but 
the Caroli-deGennes-Matricon mode. 3 ) We could ex- 
pect from the different character in quasiparticles that 
the flux flow conductivity rjf 4 ' is different between dirty 
and clean superconductors. This expectation is not ful- 
filled in isotropic s-wave superconductors; <7f in clean 
isotropic s-wave superconductors 5,6 -' turns out to be dif- 
ferent from that 7 ) in dirty superconductors only by a 
factor of In (Aqo/T) at temperature T, where Aoo is the 
modulus of the pair-potential in the spatially uniform 
state. 

However, this is not the whole story. Volovik 8 ) pointed 
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out that a single impurity in vortex cores of chiral p- 
wave superconductors 9 ^ does not change the spectrum. 
At low temperatures, the impurity scattering is expected 
to be the main process of relaxation. Further, the relax- 
ation time governs the conductivity in vortex states. 10, n ) 
Therefore, if the impurity scattering rate of quasiparti- 
cles in chiral p-wave vortex is completely different from 
those in isotropic s-wave vortex, we expect different be- 
haviors of flux flow conductivity between chiral p-wave 
and isotropic s-wave superconductors. 

The chiral p-wave superconducting state is specified 
by the d-vector 12, 13 ) d — z(p x ±ip y ) 9 ^ and expected 
to be realized in Sr? Ru0 4 . 14) The chiral p-wave state 
has two-fold degeneracy. Each phase is specified by the 
"chirality" , which we define here as the angular momen- 
tum of relative motion of Cooper-pairs; the phase with 
d = z [p x + ip y ] (z [p x — ip y ]) has +1 (—1) chirality. In 
the absence of external magnetic field, two homogeneous 
phases with definite chirality are energetically equal 9 ) 
and hence it depends on the history what phase is real- 
ized among the following possibilities: spatially uniform 
phase with positive or negative chirality or the domain 
structure consisting of phases with different chiralities. 

The magnetic field lifts the degeneracy of the two ho- 
mogeneous phases. 15 -* If the external field induces vor- 
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tices with positive vorticity, the phase with negative chi- 
rality becomes stable and that with positive chirality be- 
comes metastable. When the system is axisymmetric, 
the center-of-mass motion of Cooper-pairs acquires an- 
gular momentum (vorticity) with respect to the vortex 
center. Total angular momentum is given by the sum 
of vorticity and chirality. In terms of the total angular 
momentum around vortex center, the phase with L z = 2 
is metastable and the phase with L z = is stable. These 
two phases are schematically described in Fig. 1. 

Although the degeneracy is lifted in vortex state, 
Ginzburg-Landau theory gives small difference of energy 
between the two states. 15 ) Therefore we expect that the 
phase L z = 2 is relevant to experiments through the do- 
main phase. We expect that the phase with L z = is 
realized in field-cooling experiments and domain struc- 
ture consisting of the phase with L z = and L z = 2 
is realized in zcro-ficld-cooling experiments . We con- 
sider, in the following part, single vortex L z = and 
L z = 2. We also consider isotropic s-wave vortex, which 
has L z = 1, as a reference. 




Lz=2 Lz=0 



Fig. 1. Schematic description of chiral p-wave vortex. Cooper- 
pairs of chiral p-wave vortex states have ±1 angular momentum 
( chirality ) in the internal motion. In the presence of vortex, 
the center-of-mass motion of Cooper-pairs has the angular mo- 
mentum (vorticity) around a vortex center. The total angular 
momentum is given by the sum of chirality and vorticity 

By examining the impurity scattering rates, or the in- 
verse of impurity relaxation time, we can find the ev- 
idence that quasiparticles in vortex cores are different 
from that in normal states. The energy dependence of T 
(normalized by the normal state impurity scattering rate 
T n ) of quasiparticles (Andreev bound states) localized in 
vortex cores is given by 

f ci In (A^/E) , for L z = 1 (s-wave) 11 ) 
r/r n = I c 2 , for L z = 2 (chiral p-wave) 16 ) 

[ 0, for L z = (chiral p-wave) 16 ) 

(1.1) 

within the non-selfconsistent Born approximation 11 ) or 
the self-consistent Born approximation with a restric- 
tion on the energy range. 16 ' 17 ) Here ci and c 2 are 
independent of energy. These two are of the order 



of 0(1/ In (Aqo/T)) when the vortex shrinking (the 
Kramer-Pesch effect 18 )) occurs; otherwise they are of 
the order of unity. The phase of the pair-potential de- 
pends on the pairing symmetry. Further the sensitivity 
on the phase of the pair-potential is characteristic of the 
Andreev bound states. Therefore, the impurity scatter- 
ing rate of the Andreev bound states depends on the 
pairing symmetry. From (1.1), we can see that in s- 
wave vortex, the bound states with smaller energy are 
scattered more strongly by impurities. In chiral p-wave 
vortex with L z = 2, the bound states with various en- 
ergy are scattered equally. In chiral p-wave vortex with 
L z = 0, the bound states are not subject to the impu- 
rity scattering; this point is worthy of attention. The 
cancellation of angular momentum for a chiral p-wave 
vortex makes the physics in vortex cores similar to that 
in s-wave superconductors in zero magnetic field. We 
know that non-magnetic impurities do not affect super- 
conducting properties in the s-wave superconductors in 
the absence of magnetic field (the Anderson theorem 19 ) ) . 
In ref. 20, the authors have interpreted the vanishing T 
in chiral p-wave vortices with L z = as a consequence 
of novel applicability of the Anderson theorem in those 
vortex cores. 

The difference of the energy dependence of T (1.1) 
leads to different temperature dependence of flux flow 
conductivity for the three pairing symmetries; with the 
basis of (1.1) and consideration of the vortex shrinking 
effect (the Kramer-Pesch effect), the flux flow conductiv- 
ity CTf normalized by the normal state conductivity a n is 
expected to be 4, 16 ) 

Of/cr n 

r ~ In (Aoo/T) H c2 /B, for L z = 1 (s-wave) 6 ) 

< - [In (Aoo/T)] 2 H c2 /B, for L z = 2 (chiral p-wave^) 

{ > C(A oc /T)i? c2 / J B, for L z = (chiral p-wave) 16 ) 

at temperature T and under the magnetic field B for 
clean superconductors with the upper critical field H c2 . 
Thus, it is of great importance to confirm the result (1.1) 
on the energy dependence of impurity scattering rate in 
order to discuss the flux flow (1.2). 

In this paper, we examine the results (1.1) by numer- 
ical calculations. The scheme we adopt is the quasi- 
classical theory of superconductivity 2 1_24 ) and the self- 
consistent Born approximation. Before going into the de- 
scription of model, methods and other details, we present 
main result in this paper. Figure 2 shows the impurity 
scattering rate as a function of energy of Andreev bound 
states. The impurity scattering rate within an s-wave 
vortex with r n / Aqo = 0.01 shown by solid circles re- 
veals the logarithmic energy dependence. For chiral p- 
wave vortex with L z = 2 and r n /Aoo = 0.01 (open cir- 
cles), the impurity scattering rate is less dependent on 
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energy and particularly becomes independent of energy 
for E/Aoo < 0.1. For chiral p-wave vortex with L z = 
and Tn/Aoo = 0.1 (solid squares), the impurity scatter- 
ing rate is extremely small. From these observations, we 
confirm the previous result 11 ' 16 - ) (1.1). Accordingly, the 
result on the flux flow conductivity (1.2) becomes more 
convincing. 




Fig. 2. The impurity scattering rate T as a function of energy of 
the Andreev bound states. Solid circles represent the results for 
s-wave vortex with r n /Aoo = 0.01, open circles for chiral p-wave 
vortex with L z = 2 and r n /Aoo = 0.01, and open squares for 
chiral p-wavc vortex with L z = and Tn/Aoo = 0.1 



The rest of this paper is organized as follows. In 
the next section, we first summarize the formulation of 
quasiclassical theory of superconductivity 21-24 -* and self- 
consistent Born approximation, next explain the way 
how we can investigate the bound state in the presence of 
impurities within the green function formalism and third 
give a brief description on numerics. In section III, we 
present numerical results on the angular resolved local 
density of states for the three pairing symmetries. From 
those results, we can obtain the impurity scattering rate 
shown in Fig. 2. In section IV, we discuss the signifi- 
cance of our findings, in comparison with earlier studies. 
In section V, we summarize our conclusions. 

§2. Model and Method 

2. 1 quasiclassical theory of superconductivity 

We study the impurity effects in vortex cores in two- 
dimensional superconductors with s-wave and chiral p- 
wave pairing symmetries. For simplicity, we consider the 
systems in the type II limit with the circular symmetric 
Fermi surface. In the quasiclassical theory of supercon- 
ductivity, s-wave and chiral p-wave superconductors can 
be studied on equal- footing. In what follows, we adhere 
to the notations in rcf. 16. What we calculate and dis- 
cuss in this paper is the retarded part of the quasiclassical 



green function in the equilibrium case 
g(e,r,i>) = g(e,r,a) = 



9 f 

f -9 



(2.1) 



which is a 2 x 2 matrix in particle-hole space and is a 
function of frequency e, the direction p — (cos a, sin a) 
of momentum p — ppp and the point r = r(cos (f>, sin 4>) 
in real space. The quasiclassical green function satisfies 
the normalization condition 



-7r 2 i (1; the 2x2 unit matrix), 



(2.2) 



and follows the equation of motion (the Eilenberger equa- 
tion) 21 ) 



iv ■ Vg = (e + i5)f 3 - A - Y,,g 



(2.3) 



Here v = vp is the Fermi velocity, S is a positive in- 
finitesimal and f 3 is a Pauli matrix 



1 

-1 



The matrix A = A(r,p) in eq. (2.3) is given by 



A(r,p) 
-A*(r,p) 



(2.4) 



(2.5) 



in terms of the pair-potential A(r,p). The p (or equiva- 
lently a) dependence of the pair-potential is determined 
by the pairing symmetry The self-energy £ = S(e, r) 
in the right hand side of (2.3) is given by that in the 
self-consistent Born approximation 



T n (g(e,r,a)), 



(2.6) 



where T n is the impurity scattering rate in the normal 
state and (■••) = / • • • da/ (2ir) denotes the average over 
the Fermi surface. 

We assume that the pair-potential A(r,p) has the fol- 
lowing form: 



A(r,p) = A (r)e i (^ Q ' +ii " a . 



(2.7) 



Here Ao(r) is a monotonically increasing function sat- 
isfying the following conditions: Ao(r = 0) = and 
lim^oo Ao(r) = Aoo. The expression (2.7) describes iso- 
lated vortices with vorticity +1 at r = 0. Under a given 
form of the pair-potential, g and £ are self-consistently 
determined in numerical calculations. For a given p or 
a, a rotated frame in real space 



r = sp + bz x p 



(2.8) 



with s = cos (0 — a) and b — r sin ((f) — a) is suitable to 
study eq. (2.3). The left hand side of eq. (2.3) can be 
rewritten as —ivdg/ds and the equation yields a one- 
dimensional problem on the line (the quasiclassical tra- 
jectory) with constant b. Now s and b turn out to be, 
respectively, the coordinate along the quasiclassical tra- 
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jectory and the impact parameter with respect to the 
vortex center. 



and E p (b) is given by 




Fig. 3. Quasiclassical trajectories, on which the Eilenberger 
equation is to be solved, are shown for a given direction a of 
momentum of quasiparticlcs. Solid circles represent the points 
where wc calculate the angular-resolved local density of states. 



2.2 Andreev bound states and angular resolved local 
density of states 
For energy much lower than the modulus of the pair- 
potential Aqo in the bulk, most spectral weight of g is 
exhausted by the contribution with the impact param- 
eter b whose modulus is much smaller than the coher- 
ence length £ = v/^Aoo). In the absence of impurities 
(and hence the self-energy in eq. (2.3)), the quasiclassical 
green function for |e| -C Aoc and \b\ <C £o> is expected to 
have the form of 4 ' 16 - 18 < 25 ) 

Trvexp [— u(s)] Mia) , , , ,„ „. 

g < C ' r ' a > = C(e-^ p (6) + iV +(reSUlarPart) ' ^ 
as a function of e, s, b and a. Here the function u(s) in 



the numerator is defined by 



u(s) = -/*'' ds'Ao(s') (2.10) 
v Jo 

and the matrix M(a) is defined by 

/ 1 -ic iLi " \ 

M(a) = >\ . (2.11) 



The constant C in the denominator denotes the normal- 
ization factor 

/>oo 

C= / dsexp[-u(s)] - Co (2.12) 
Jo 



b 

C 



' A («) _ 
as e 



X( S ) 



(2.13) 



The subscript p stands for pure superconductors. The 
second term in the right hand side of (2.9) represents neg- 
ligible contributions from the scattering states of quasi- 
particles. 

From the above results, we see that N(e, a, s, b)/No — 
-Img(e,r,ct), which is "the angular-resolved local den- 
sity of states" 25 ~ 27 ) normalized by the density of states 
No in the normal state, has the following expression: 

N{e, a, s, b)/N = ^vC^e'^ 8{e - E p (b)). (2.14) 

This expression (2.14) has a sharp peak along a tra- 
jectory shown in Fig. 3; on the trajectory, the relation 
e = E p (b) is satisfied for a given a. Numerical results on 
the angular-resolved local density of states were firstly 
presented by Klein 26 ) for vortex lattice of an s-wave su- 
perconductor without impurities. For isolated s-wave 
vortex with impurities, the angular resoloved local den- 
sity of states were presented in ref. 25. 

In the presence of sufficiently dilute impurities, the 
expression (2.9) is expected to be modified as 16 ' 17 ) 



7ruexp [— u(s)] M(a) 



+ (regular part), 



C (e - E(b) + iT) 
and accordingly the expression (2.14) turns into 

i e -«(s)p 



N(e, a, s, b)/No = 



(2.15) 



(2.16) 



[(e-E(b)r+^Y 

where E(b) and T are certain functions of b and are in- 
dependent of s for a given a. These results are sug- 
gestive of the way to investigate the energy of Andreev 
bound states and the impurity scattering rate in impure 
but clean superconductors; if N(e,a,s,b)/N has a sin- 
gle peak as a function of e, the peak position gives the 
energy E(b) of Andreev bound states for a given impact 
parameter b. The width T of the single peak, on the other 
hand, gives the impurity scattering rate for the Andreev 
bound states. In the next section, we will present nu- 
merical results on N(e,a,s,b)/N as functions of e for 
a certain fixed value of a, s(= 0) and various values 
of b, or cquivalently for the points shown by solid cir- 
cles in Fig. 3. Furthermore, owing to the axisymmetry, 
N(e, a, s — 0,b)/N o is independent of a. We can thus 
set a = without loss of generality 

2.3 Numerical details 

We take the profile of the pair-potential as 



A (r) = Aootanhtr/fo)- 



(2.17) 



The Eilenberger equation (2.3) is solved numerically by 
the fourth-order Rungc-Kutta method after transform- 
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ing to the Riccati formulation 25 ' 28 ~ 32 ) of (2.3). The cut- 
off radius, from which the Runge-Kutta program starts, 
is taken as 3£o- In the Runge-Kutta method, the choice 
of the increment step being 3 x 10~ 3 £o along the qua- 
siclassical trajectory yields the bounds of relative error 
10~ 4 . The a-intcgration on the 2D Fermi surface is re- 
placed by the sum over the points ctj = 2irj/N a with 
j 6 [1, N a \. The choice of N a = 3000 for T n = O.lAoo is 
sufficient to obtain the upper bound of the error of the 
self-energy ST, 



Max(|js n |,|<ys 12 |,|js 21 |) < io- 3 r„ 



(2.18) 



In our calculation, the self-consistent condition on the 
self-energy is satisfied within the accuracy (2.18). 

§3. Results 

3.1 S-wave vortex 
Figure 4 shows e dependence of the angular-resolved 





Fig. 4. The angular resolved local density of states 

N(e,a = 0, s = 0,b)/No of s-wave superconductors for fc/£o = 
0.01 and 0.05n with n = 1,2, ■ • • for (a) T n = O.OlAoo and (b) 
r n = O.lAoo. For both cases, the smaller the impact parameter 
is, the width of the peak becomes broader. 



local density of states of the quasiparticles with momen- 
tum direction a = and at the positions (shown by 



solid circles in Fig. 3) of s = and fe/£o = 0.01 and 
0.05n with n = 1, 2, • • • (from left to right). The upper 
panel (a) represents the results for T n = O.OIAqo and 
the lower panel (b) for T n = O.lAoo- Each curve for a 
given impact parameter b has a single peak; This fact 
shows that the picture of the bound states still works 
even in the presence of impurities. First we discuss 
the peak position. The smaller the impact parameter 
b is, the peak position shifts to the lower energy. Fig- 
ure 5 shows that the energy of bound states for a given 
impact parameter for T n = O.OIAqo (solid circle) and 
r n = O.lAoo (open circle). The line represents analyti- 
cal result E p (b) — 0. 852557 Aoo^/Co, which is obtained by 
substituting (2.17) into (2.13) for pure superconductors. 
The bound state energy is nearly linear in the impact 
parameter. This bound state energy for T n = O.OlAoo is 
almost same as that of pure case. The bound state en- 
ergy for r n = O.IAqo is enhanced compared to the pure 
case, owing to the self-energy effect. 

Next we discuss the width of each curve in fig. 4. In 
the figures, the smaller the impact parameter is, the peak 
becomes broader. This means that the Andreev bound 
states with smaller energy are scattered more strongly 
by impurities. This observation agrees qualitatively with 
the results of non-self consistent calculations mentioned 
in the Introduction. For T n = O.OlAoo, each curve fits 
well to the form of the Lorentzian 1/ (e — E) 2 + V 2 up 
to a constant factor. From this fitting, we obtain the 
impurity scattering rate as a function of bound state 
energy, which has been already shown in Fig. 2. For 
T n = O.lAoo, on the other hand, each curve does not 
necessarily fit to the Lorentzian; a comparison of the im- 
purity scattering rate calculated from the curvature near 
the peak position and that from the full width of the 
half-height shows appreciable difference, particularly for 
small impact parameters. Therefore, we leave the dis- 
cussion for T n = O.lAoo in a qualitative level. 



3. 2 Chiral p-wave vortices with L z = 2 

Figure 6 shows e dependence of the angular-resolved 
local density of states for chiral p-wave vortices with 
L z = 2 for (a) T n = O.OlAoo and (b) T n = O.lAoo- 
The parameters such as a, s and b are taken as the same 
as that in the s-wave case. In figures, more right curve 
corresponds to larger impact parameter. Single peak for 
each impact parameter suggests the existence of a bound 
state. The bound state energy obtained from the peak 
position is shown in fig. 7. The bound state energy for 
r n = O.OlAoo shown by solid circles is almost same as 
that of pure superconductors shown by the solid line. 
On the other hand, the bound state energy T n = O.lAoo 
shown by open circles is suppressed compared to the pure 
case. This tendency of the renormalization effect is oppo- 
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Fig. 5. The bound state energy as a function of the impact pa- 
rameter for s-wave vortex. Solid circles represent the results for 
r n = O.OlAoo and open circles T n = O.lAoo. The line represents 
analytical result E p (b) = 0.852557Aoo b/£o, which is obtained by 
substituting (2.17) into (2.13) for pure superconductors. 

site to that of s-wave case, where the self-energy renor- 
malization yields the enhancement of the bound state 
energy. 

Now we turn to the width of the peak. The shapes 
of the curves are almost same for small impact param- 
eters. This means that the impurity scattering rate be- 
comes independent of the impact parameter b for small b 
or equivalcntly the bound state energy for small energy 
This is confirmed quantitatively; for r n /Aoo = 0.01, all 
curves fit the Lorentian. The deduced impurity scatter- 
ing rate has been shown in fig. 2. Our results confirm the 
previous results (1.1) for chiral p-wave case with L z = 2. 

3.3 Chiral p-wave vortices with L z = 

Figure 8 shows e dependence of the angular-resolved 
local density of states for chiral p-wave vortices with 
L z = for T n = O.lAoo (Numerical calculation for 
T n = O.OlAoo turns out to be prohibitive). The param- 
eters such as a, s and b are taken as the same as that in 
the s-wave case. In figures, more right curve corresponds 
to larger impact parameter. Open circles in fig. 9 rep- 
resent the bound state energy obtained from the peak 
position. The solid line shows the bound state energy 
in pure superconductors. The self-energy leads to the 
suppression of the bound state energy, contrary to the s- 
wave case but similarly to the chiral p-wave vortex with 
L z = 2. 

Now we turn to the width of the peak, which is obvi- 
ously much smaller than those for s-wave (Fig. 4b) and 
chiral p-wave with L z = 2 (Fig. 6b). The Lorentian-fit 
works well for all curves and the impurity scattering rate 
r is deduced, as shown in Fig. 2. Finite but extremely 
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Fig. 6. The angular resolved local density of states 

N(e, a = 0, s = 0, b) of chiral p-wave superconductors with L z = 
2 for 6/£ = 0.01 and 0.05n with n = 1, 2, • • • for (a) T n = 
O.OlAoo and (b) T n = O.lAoo. For both cases, the width of all 
the curves is almost same, in contrast to the s-wave. 
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Fig. 7. The bound state energy as a function of the impact pa- 
rameter for chiral p-wave vortex. Solid circles represent the re- 
sults for T n = O.OlAoo and open circles T n = O.lAoo- The line 
represents analytical result E p (b) = 0.852557Aoob/£o, which is 
obtained by substituting (2.17) into (2.13) for pure superconduc- 
tors. 



small impurity scattering rate is consistent with the pre- 
vious results (1.1) for chiral p-wave with L z = 0. 

§4. Discussion 

Here we recall earlier studies on the impurity problem 
of chiral p-wave vortex and make clear the importance of 
our findings. Volovik 8 ) studied energy spectrum of quasi- 
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Fig. 8. The angular resolved local density of states 

N(e, a = 0, s = 0, 6) of chiral p-wave superconductors with L z = 
for b/£ = 0.01 and 0.05n with n = 1, 2, • • • for T n = O.lAoo. 
Peak width is much smaller than those in s-wave vortex(Fig.4b) 
and chiral p-wave vortex with L z = 2 (Fig. 6b) with the same 
value of T n . 




Fig. 9. The bound state energy as a function of the impact pa- 
rameter for chiral p-wave vortex. Open circles represent the 
results for T n = O.lAoo. The solid line represents analytical re- 
sult E p (b) = 0.852557A 00 b/£o, which is obtained by substituting 
(2.17) into (2.13) for pure superconductors. 



particles bound to vortex cores. His method is a com- 
bination of the Andreev approximation of Bogoliubov- 
deGennes (BdG) equation and Bohr-Sommerfeld semi- 
classical quantization. He concluded that a single im- 
purity does not change the spectrum of chiral p-wave 
vortex with L z = 2 and L z = 0; it is crucial, in his 
theory, whether the chirality is odd or even. However, 
this statement is subject to an counterexample. Mat- 
sumoto and Sigrist 33 ' found in numerical calculation of 
BdG equation of chiral p-wave vortex that the single im- 
purity at the vortex center shifts the energy spectrum 
of Caroli-deGennes-Matricon mode for L z = 2 and does 



not for L z = 0. Very recently, Miyazu 34 ' studied analyt- 
ically the effect of single impurity at a generic position 
inside cores of chiral p-wave vortex and found that the 
impurity is harmless for L z = and effective for L z = 2. 

The above results are on single impurity in chiral p- 
wave vortex. However, in realistic situations there are 
many impurities. Further, there are no explicit impli- 
cations on physical quantities accessible in experiments. 
Our previous results 16 ' in (1.1) for chiral p-wave vor- 
tex have addressed these issues. The results in the 
present paper further confirm that the results (1.1) are 
correct even in the self- consistent treatment of impurity 
self-energy. Accordingly we expect that the prediction 
on the flux flow conductivity (1.2) comes true. The 
measurements of flux flow conductivity of chiral p-wave 
superconductors (e.g. S^RuO^) are highly desirable. 
The absence of impurity effects in vortex cores has an- 
other implication; the absence of impurity effects makes 
the Kramcr-Pesch effect much more accessible in exper- 
iments of chiral p-wave superconductors than any other 
superconductors. 20 ' This can be, in principle, examined 
experimentally, e.g. by Muon spin relaxation measure- 
ments. 35 ' Further, the cancellation of angular momen- 
tum in chiral p-wave vortex yields novel quantum effects 
on the pinning problem due to impurities. 36 ' This phe- 
nomenon can be also tested in the measurements of M-H 
curves in multiple cycles. 

§5. Conclusion 

We studied the impurity effects numerically on quasi- 
particles localized in vortex cores of s-wave and chiral 
p-wave superconductors under the scheme of quasiclassi- 
cal theory of superconductivity and self-consistent Born 
approximation. From angular-resolved local density of 
states, we calculated the bound state energy and impu- 
rity scattering rates. In s-wave case, the bound states 
with smaller energy are scattered more strongly by im- 
purities. On the other hand, in the case of the chiral 
p-wave vortex with L z = 2, the impurity scattering rate 
is independent of energy of the bound state. An intrigu- 
ing thing occurs in the case of chiral p-wave vortex with 
L z = 0; the impurity scattering rate is almost vanishing. 
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